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Abstract. Probabilities of collisions of migrating small bodies and dust particles produced by 
these bodies with planets were studied. Various Jupiter-family comets, Halley-type comets, long- 
period comets, trans-Neptunian objects, and asteroids were considered. The total probability of 
collisions of any considered body or particle with all planets did not exceed 0.2. The amount of 
water delivered from outside of Jupiter's orbit to the Earth during the formation of the giant 
planets could exceed the amount of water in Earth's oceans. The ratio of the mass of water 
^sO , delivered to a planet by Jupiter-family comets or Halley-type comets to the mass of the planet 
can be greater for Mars, Venus, and Mercury, than that for Earth. 
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1. Model of migration of small bodies and dust particles 



In the present paper we discuss the probabilities of collisions of migrating small bodies 
and dust particles produced by these bodies with planets. Ipatov & Mather (2003, 2004a- 
b, 2006, 2007), Ipatov et al. (2004), and several other authors cited in the above papers 
studied the probabilities with Earth, Venus, and Mars. In our above papers (which can 
be found on astro-ph and http://faculty.cua.edu/ipatov/list-publications.htm), the prob- 
abilities of the collisions were calculated based on the time variations of orbits of bodies 
and particles during their dynamical lifetimes (until all bodies or particles reached 2000 
AU from the Sun or collided with the Sun). Using the same time variations, below we 
consider also the probabilities of collisions with the giant planets and Mercury. 
• The orbital evolution of >30,000 bodies with initial orbits close to those of Jupiter- 



family comets (JFCs), Halley-type comets, long-period comets, trans- Neptunian objects, 
and asteroids in the resonances 3/1 and 5/2 with Jupiter, and also of >20,000 dust 
Q\ ' particles produced by these small bodies was integrated during their dynamical lifetimes. 
■ We considered the gravitational influence of planets, but omitted the influence of Mercury 
(exclusive for Comet 2P/Encke). In about a half of calculations of migration of bodies, we 
used the method by Bulirsh-Stoer (1966) (BULSTO code), and in other runs we used a 
symplectic method (RMVS3 code). The integration package of Levison & Duncan (1994) 
was used. For dust particles, only the BULSTO code was used, and the gravitational 
influence of all planets, the Poynting-Robertson drag, radiation pressure, and solar wind 
drag were taken into account. The ratio (3 between the radiation pressure force and the 
gravitational force varied from ^0.0004 to 0.4. For silicates, such values of (3 correspond 
to particle diameters d between ^1000 and 1 microns; d is proportional to 1//3. 

In our calculations, planets were considered as material points, so literal collisions did 
not occur. However, using the algorithm suggested by Ipatov (1988) with the correction 
that takes into account a different velocity at different parts of the orbit (Ipatov & Mather 
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2003), and based on all orbital elements sampled with a 10-500 yr step, we calculated 
the mean probability P of collisions of migrating objects with a planet. The step could 
be different for different typical dynamical lifetimes of particles. We define P as P^/N, 
where is the probability of a collision of all N objects with a planet. The probabilities 
of collisions of bodies and particles at different (3 with planets are presented in Fig. 1. 
These probabilities do not take into account the destruction of particles in collisions and 
sublimation, which can be more important for small particles. Our runs were made mainly 
for direct modelling of collisions with the Sun, but Ipatov & Mather (2007) obtained that 
the mean probabilities of collisions of considered bodies with planets, lifetimes of the 
bodies that spent millions of years in Earth-crossing orbits, and other obtained results 
were practically the same if we consider that bodies disappear when perihelion distance 
becomes less than the radius of the Sun or even several such radii. 

2. Probabilities of collisions of migrating small bodies and dust 
particles with planets and the Sun 

The probability Pe of a collision of a JFC with the Earth exceeded 4 • 10~ 6 if initial 
orbits of bodies were close to those of several tens of JFCs, even excluding a few bodies 
for which the probability of a collision of one body with the Earth could be greater 
than the sum of probabilities for thousands of other bodies. The Bulirsh-Stoer method 
of integration and the symplectic method gave similar results. The ratios of probabilities 
of collisions of JFCs with Venus, Mars, and Mercury to the mass of a planet usually 
were not smaller than those for Earth. For most considered bodies, the probabilities 
Pue of collisions of most bodies with Mercury (exclusive for Comet 2P /Encke, for which 
Pme^Pe) were smaller by an order of magnitude than those with Earth or Venus. 

For dust particles produced by comets and asteroids, Pe was found to have a maximum 
(-0.001-0.02) at 0.002 < (3 < 0.01, i.e., at d~100 fim (this value of d is in accordance 
with observational data) . These maximum values of Pe were usually (exclusive for Comet 
2P /Encke) greater at least by an order of magnitude than the values for parent comets. 
Probabilities of collisions of most considered particles with Venus did not differ much 
from those with Earth, and those with Mars were about an order of magnitude smaller. 
For particles produced by Halley-type comets, P was greater for Mercury than for Mars. 

Using Pe = 4 • 10~ 6 and assuming that the total mass of planetesimals that ever 
crossed Jupiter's orbit was about lOOm^ (Ipatov 1987, 1993), where vtle is the mass 
of the Earth, Ipatov & Mather (2003, 2004a-b, 2007) concluded that the total mass of 
water delivered from the feeding zone of the giant planets to the Earth could be about 
the total mass of water in Earth's oceans. (Similar conclusion was made by Ipatov (2001) 
based on other calculations.) We supposed that the fraction of water in planetesimals 
equaled 0.5. The ratio of the mass of water delivered to a planet by Jupiter-family comets 
and Halley-type comets to the mass of the planet can be greater for Mars, Venus, and 
Mercury, than that for Earth. This larger mass fraction would result in relatively large 
ancient oceans on Mars and Venus. The larger value of P for Earth we have calculated 
compared to those argued by Morbidclli et al. (2000) (Pe — (1 — 3) • 10~ 6 ) and Levison et 
al. (2001) (Pe — 4- 10~ 7 ) is caused by the fact that Levison et al. (2001) did not take into 
account the gravitational influence of the terrestrial planets, and Morbidclli et al. (2000) 
considered low-eccentric initial orbits beyond Jupiter's orbit. Besides, we considered a 
larger number of bodies. The detailed discussion on delivery of water and the comparison 
of our results with the results by other authors were made by Ipatov & Mather (2007). 
Marov & Ipatov (2005) discussed the delivery of volatiles to the terrestrial planets. 

At the present time, most authors (e.g., Lunine et al. 2003, Morbidelli et al. 2000, and 
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Petit et al. 2001) consider that the outer asteroid belt was the main source of the delivery 
of water to the terrestrial planets. Drake & Campins (2006) noted that the key argu- 
ment against an asteroidal source of Earth's water is that the O's isotopic composition 
of Earth's primitive upper mantle matches that of anhydrous ordinary chondrites, not 
hydrous carbonaceous chondrites. To the discussion of the deuterium/hydrogen paradox 
of the Earth's oceans presented by Ipatov & Mather (2007), we can add that Genda & 
Ikoma (2008) showed that D/H in the Earth's ocean increased by a factor of 2-9. 

Probabilities of collisions of JFCs with Saturn typically were smaller by an order of 
magnitude than those with Jupiter, and collision probabilities with Uranus and Neptune 
typically were smaller by three orders of magnitude than those with Jupiter. As only a 
small fraction of comets collided with all planets during dynamical lifetimes of comets, 
the orbital evolution of comets for the considered model of material points was close to 
that for the model when comets collided with a planet are removed from integrations. 

Probabilities of collisions of considered particles and bodies with Jupiter during their 
dynamical lifetimes are smaller than 0.1. They can reach 0.01-0.1 for bodies and particles 
initially moved beyond Jupiter's orbit or in Encke-type orbits. For bodies and particles 
initially moved inside Jupiter's orbit, the probabilities are usually smaller than the above 
range and can be zero. Probabilities of collisions of migrating particles (exclusive for 
trans-Neptunian particles) with other giant planets were usually smaller than those with 
Jupiter. The total probability of collisions of any considered body or particle with all 
planets did not exceed 0.2. 

Collisions of planetesimals with a star can cause variations in observed brightness and 
spectrum of the star. In our calculations, the fraction Ps un of comets collided with the 
Sun during their dynamical lifetimes was about a few percent. For most JFCs, dynamical 
lifetimes are less than 10 Myr, and on average P5„„^0.02. For dust particles, Ps un 
depends on (3 and can be considerably greater than for their parent bodies. For example, 
for Comet lOP/Tempel 2, Ps„„«0.01, and almost all particles produced by this comet 
collide with the Sun (Ipatov & Mather 2006). 
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Caption to the Figure: 

The probability of collisions of dust particles and bodies (during their dynamical life- 
times) with Mercury (subfigure Me), Venus (V), Earth (E), Mars (Ma), Jupiter (J), 
Saturn (S), Uranus (U), and Neptune (N) versus (3 (the ratio of the radiation to gravita- 
tional forces) for particles launched from asteroids (ast), trans-Neptunian objects (tno), 
Comet 2P/Encke at perihelion (2P per), Comet 2P/Encke at aphelion (2P aph), Comet 
2P/Enckc in the middle between perihelion and aphelion (2P m), Comet lOP/Tempel 2 
(10P), Comet 39P/Oterma (39P), long-period comets (Ip) at eccentricity e=0.995 and 
perihelion distance q=0.9 AU, and Halley-type comets (ht) at e=0.975 and g=0.5 AU 
(for Ip and ht runs, initial inclinations were from to 180°, and particles were launched 
near perihelia). If there are two points for the same /3, then a plot is drawn via their 
mean value. Probabilities presented at /3^10~ 5 arc for small bodies (/3=0). Probabili- 
ties presented only for bodies were calculated for initial orbits close to orbits of Comets 
9P/Tempel 1 (9P), 22P/Kopff (22P), 28P/Neujmin (28P), 44P/Reinmuth 2 (44P), and 
test asteroids from resonances 3:1 and 5:2 with Jupiter at e=0.15 and «=10° ('ast 3:1' 
and 'ast 5:2'). For series nl and n2, initial orbits of bodies were close to 10-20 different 
Jupiter- family comets (Ipatov & Mather 2004b). 
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